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ABSTRACT 



We present arguments based on the measured abundances in individual damped Lyman a systems (DLAs) and sub-damped Lyman 
a systems (sub-DLAs), and also the average abundances inferred in large samples of QSO absorption line systems, to suggest that 
the amount of dust in intervening QSO absorbers is small and is not responsible for missing many QSOs in magnitude limited QSO 
surveys. While we can not totally rule out a bimodal dust distribution with a population of very dusty, metal rich, absorbers which 
push the background QSOs below the observational threshold of current optical spectroscopic studies, based upon the current samples 
it appears that the metallicity in QSO absorbers decreases with increase in H I column densities beyond 10" cm"^. Thus the sub-DLA 
population is more metal rich than the DLAs, a trend which may possibly extend to the non-damped Lyman limit systems (NDLLS). 
Based on the recently discovered mass-metallicity relation for galaxies, we suggest that most sub-DLAs and possibly NDLLS, are 
associated with massive spiral/elliptical galaxies while most DLAs are associated with low mass galaxies. The sub-DLA galaxies will 
then contribute a larger fraction of total mass (stellar and ISM) and therefore metals, to the cosmic budget, specially at low redshifts, 
as compared to the DLAs. 

Key words. Quasars: absorption lines — ISM: abundances, dust, extinction 



1. Introduction 

Damped Lyman a systems (DLAs), having column density of 
^ ' H I (Nh i) > 2 X 10^" cm"^ have been suggested to be the high- 
redshift analogs of disks of nearby luminous galaxies. In spite 
of having a high amount of neutral hydrogen these absorbers 
have very little H2 (Ledoux et al. 2003), which raises questions 
about their being associated with high star formation activity. A 
large fraction of DLAs do indeed appear to have low star for- 
mation rates based on deep emission-line imaging searches (e.g. 
KuUcarni et al. 2006a and references therein). Metallicity mea- 
surements, using intermediate to high resolution observations, 
are available for > 100 DLAs. These indicate that the metallicity 
evolution of DLAs is weak; most of the DLAs, even at Zabs ~ 
are found to be metal poor (Kulkarni et al. 2005). It is thus pos- 
sible that the DLAs may not trace the bulk of star formation in 
the Universe and therefore may not be among the leading metal 
carriers. 

Compilation of observed metallicities of individual DLA 
systems shows a trend of decreasing abundance with increas- 
ing H 1 column density (Boisse et al. 1998; Akerman et al. 2005; 
Meiring et al. 2006a). The trend possibly exists to lower H 1 col- 
umn densities covering the sub-DLAs, having Nh i between 10'^ 
and 2 x 10^*^ cm"^ (Peroux et al. 2003). A similar conclusion was 
drawn by York et al. (2006, hereafter, Y06) based on the average 
abundances inferred from the composite spectra of sub-samples 
drawn from a sample of 809, intervening Mg 11 QSO absorp- 
tion line systems (rest equivalent width of /12796 > 0.3 A and 
1 < Zabs < 1.86) compiled from the Sloan Digital Sky Survey 



(SDSS) Data Release 1 (DRl). The sub-sample with highest es- 
timated (near solar) Zn metallicity, was found to have a low aver- 
age infeiTed Nh i of ~ 10^*^ cm"^. Several arguments were given 
in Y06 to suggest that the DLAs may not be the analogs of the 
local bright galaxies, rather, the sub-DLAs and possibly the non- 
DLA Lyman limit systems (NDLLS) with Nh i between 10'^ and 
10'^ cm~^ may be associated with large galaxies. Here we con- 
sider this scenario further, using the currently available data on 
individual absorbers and the results of Y06. We consider vari- 
ous selection effects at work and discuss the implications of the 
recently established mass-metallicity relation for the nature of 
DLAs and sub-DLAs. 

In section 2 we study the dependence of the Zn abundance 
on the column density of H I, for the observations of individual 
systems taken from the literature and for the average values of 
these quantities obtained by Y06 for large samples. We present 
other relevant observational information and discuss the obscu- 
ration threshold and other selection effects to argue for the reality 
of the observed dependence. The implication of this dependence 
for the nature of the absorbing galaxies are presented in section 
3 and conclusions are presented in section 4. 

2. Dependence of metallicity on the neutral 
hydrogen column density 

2.1. Observational data 

The number of metallicity measurements of individual sub- 
DLAs has grown considerably since the compilation of Peroux 
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et al. (2003). We have compiled a sample (the literature sam- 
ple, hereafter TLS) of all the measurements of Zn abundance 
of DLAs and sub-DLAs from the literature. This sample con- 
sists of 119 DLAs with 0.1 < Zabs < 3.9 and 30 sub-DLAs 
with 0.6 < Zabs < 3.2. The DLA Zn sample is based on our 
recent HST, MMT, and VLT data (Khare et al. 2004; Kulkarni et 
al. 2005; Meiring et al. 2006a,2006b; Peroux et al. 2006b) and 
those from the literature (data compiled in Kulkarni et al. 2005, 
and more recent data from Rao et al. 2005; Akerman et al. 2005; 
Ledoux et al. 2006). For the sub-DLAs, we compile the data 
from Lu et al. (1995, 1996); Pettini et al. (1994, 1999, 2000); 
Kulkai-ni et al. (1999); Ellison & Lopez (2001); Srianand & 
Petitjean (2001); Dessauges-Zavadsky et al. (2003); Khare et al. 
(2004); Peroux et al. (2006a); Ledoux et al. (2006); and Meiring 
et al. (2006a,2006b). The DLA sample contains 68 detections 
and 51 limits, while the sub-DLA sample contains 13 detections 
and 17 limits. Most of these limits are at high redshifts. The full 
sample used here is given in Kulkarni et al (2006b). 

In Fig. 1, we have plotted [Zn/H], as a function of Nh i for 
TLS. The trend of decreasing abundance with increasing Nh i is 
clear The Spearman rank correlation test gives the probability 
of chance correlation to be 1.6 xlO"". It can be seen from the 
figure that the trend is similar for systems with Zabs < 1 -5 and 

Zabs > 1-5. 

We note that we have only used Zn abundances in Fig.l. 
There have been studies which use the Si and S abundances in 
systems where Zn abundances are not available due to observa- 
tional limitations. We, however, refrain from using these as (i) 
given the intrinsically higher abundance and higher strengths of 
the detectable absorption lines of these elements, line saturation 
effects may be important, leading to an underestimate of their 
abundances; (ii) Si may be depleted on the dust grains (up to -1.4 
dex in cold ISM); (iii) the lines of Zn II are detectable to lower 
redshifts; and (iv) lines of S are often in the Lyman a forest. 

As noted in the last section, Y06 inferred the average abun- 
dances of various elements in the composite spectra of several 
sub-samples of SDSS DRl QSOs with strong Mg II absorbers. 
The sub-samples were chosen on the basis of various properties 
of the absorption line systems and QSOs. The column densi- 
ties of several species, including Zn II, were estimated from in- 
verse variance-weighted, arithmetic mean, normalized spectra of 
individual QSOs in the samples. The geometric mean compos- 
ite spectra of the same sub-samples were compared with simi- 
lar spectra of matching (in emission redshifts and i magnitudes) 
samples of QSOs without absorption lines in their spectra, to de- 
termine the absorber rest frame extinction, E{B - V). The col- 
umn densities of neutral hydrogen were estimated by assum- 
ing the Small Magellanic cloud (SMC) dust-to-gas ratio, as the 
SMC extinction law was found to fit the composite extinction 
curves well. In Fig. 1 we have plotted the results for various 
sub-samples of Y06. Note that we have used all the sub-samples 
from Table 1 of Y06 for which the relevant data were available, 
some of which were not listed in their Table A4. Details of their 
sub-samples which are used here are given in Table 1 . The trend 
of metallicity dependence on Nh i is similar to that in TLS, the 
Spearman rank correlation test for Y06 points gives the proba- 
bility of chance correlation to be 0.098. We have included the 
main sample (#1) of Y06 in our analysis. As the rest of the sub- 
samples are drawn from this, its inclusion may cause some bias. 
We have verified that removing this sample has very little effect 
on the correlation noted above. 

Data for individual DLAs show a trend of increasing deple- 
tion of Cr with respect to Zn, and thus, higher dust-to-metal ratio, 
for higher metallicity (Ledoux et al. 2003; Akerman et al. 2005; 




J I I I I I I I I I I I I I I I I I I I I I I I I I i_ 

19.5 20 20.5 21 21.5 22 
Log (N,,) 

Fig. 1. Metalhcity ([Zn/H]) vs. log(NH i). Solid (Red) circles rep- 
resent average values for large samples as obtained by Y06. 
Open (Black) circles and (blue) stars are for Zabs < 1-5 and 
Zabs > 1-5, respectively, for TLS. One sigma error bars are 
shown. Also shown are the best fit lines (long dashed (black) 
line and short dashed (red) line for TLS and Y06 respectively) 
obtained by ignoring the limits. The top solid line (blue) repre- 
sents the empirical obscuration bias (Nznii = 1.4 x lO'-' cm"^). 



Meiring et al. 2006a). We have confirmed that a similar trend is 
shown by the results of Y06. Other distributions, e.g. [Cr/Zn] vs. 
log(Nzn ii) and [Cr/Zn] vs. log(NH i) are also similar for Y06 and 
TLS, showing that the results of Y06, which are average values 
for large samples, are consistent with the observations of indi- 
vidual systems. 

We note that the Y06 results are based on the assumption of a 
constant dust-to-gas ratio. This ratio may depend on metallicity 
in view of the anticorrelation between [Cr/Zn] and [Zn/H] (and 
in view of Fig. 2 below). We have tried to estimate the effect of 
this on the relation between [Zn/H] and Nh i, obtained by Y06, as 
follows. We assumed that the abundance of Zn, (Xzn) is propor- 
tional to N^ J. Thus E(B - V) which is proportional to Xzn Nh i 

( 1 O+ty) 1.0 

will be proportional to Nj^ j , giving Nhi E{B - y)(i o+o) 
Note that the Y06 results assume a = 0. For different values of 
a, we determined Nh i from the E(B - V) values, using the above 
relation. Using this value of H I column density we determined 
[Zn/H] and then the slope of the best fit line between [Zn/H] and 
Log(NHi). Note that by construction, this slope should equal a. 
We found that for positive values of a the slope remains nega- 
tive for a < 0.81 and beyond that remains < a/2 until a - 3. 
For negative values of a the slope is somewhat smaller than the 
assumed value of a and comes close to it (-0.94) for a = -0.9. 
Thus we see that assumption of a metallicity dependent dust-to- 
gas ratio, in fact, makes the decrease of Zn abundance with Nh i 
even steeper (slope — 0.9) than that (slope - -0.45) found for 
the assumption of constant-dust-to-gas ratio (Fig. 1). 

From an earlier version of this diagram (Fig. 1), with 37 
points, including limits, Boisse et al. (1998) noted the trend of 
decreasing abundance with increasing Nhi (their Fig. 19; here- 
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Table 1. Relevant data for Y06 sub-samples 
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Sample 


Selection Criterion'' 


Number 


LX)g(!Szn ll) 


LOg^JNn I) 


(oMU ) 


number 




of systems 






0.002 


24 


- i f <0.2 


698 


12.4 


20.0 


0.003 


23 


W <2 
Zab. <1.3127 


558 


12.2 


20.1 


0.006 


9 


404 


>12.4 


20.4 


0.007 


5 


1.53<WMgn <1.91 


139 


12.6 


20.5 


0.009 


11 


\f <19.12 


398 


>12.4 


20.6 


0.010 


13 


j8« <0.103 


405 


<12.3 


20.7 


0.011 


12 


i>19.12 


411 


>12.4 


20.7 


0.011 


14 


y8>0.103 


404 


12.5 


20.7 


0.012 


26 


WMgii >2.5, A(g - < 0.2 


97 


12.5 


20.7 


0.013 


1 


Full sample 


809 


12.3 


20.8 


0.014 


16 


W*^i,/WMgii >0.577 


369 


12.4 


20.8 


0.018 


6 


1.91<WMgn <2.52 


132 


12.8 


20.9 


0.019 


20 


Fe n /12374 present 


392 


12.5 


20.9 




/ 


Z.jZSVVMgii <J.\J 




LZ.y 


Zl.l 


0.032 


8 


Wwgii >2.0 


251 


12.6 


21.2 


0.034 


17 


W AV^ <1 538 


85 


12.5 


21.2 


0.034 


19 


Fe 11 /t2260 present 


58 


13.0 


21.2 


0.036 


21 


Zn II-Mg I /12026 present 


83 


13.1 


21.2 


0.058 


22 


Zn Il-Cr n i2062 present 


31 


13.3 


21.4 


0.081 


25 


A(g - >0.2 


111 


>12.7 


21.6 


0.085 


27 


WMgii >2.5,Afe-0 >0.2 


48 


>12.4 


21.6 



" For selection from the full sample (sample 1) 

* Estimated from equivalent widths in composite spectra after correction for blends 
Estimated from E(B-V) assuming a constant dust-to-gas ratio 

The difference between the actual colours of QSOs and the median colours of QSOs at that 

redshift (Richards et al. 2003) 
" Rest equivalent width of Mg II /12796 in A 
^ i magnitude of the QSO 

* Relative velocity w.r.t. the QSO 

Rest equivalent width of Fe II /!2382 in A 
' Rest equivalent width of Al II /11670 in A 
j Rest equivalent width of Mg I A2S52 in A 



after, the Boisse plot). They interpreted this as being due to 
the observational limitations in detecting weak Zn 11 lines at 
the low Nh I end and an obscuration bias (caused by dimming 
of the QSOs, by the dust in the absorbers, below the limit of 
magnitude limited surveys) towards the high Nh i end, causing 
a dearth of points in the lower left and upper right comers of 
the plot respectively. Boisse et al. (1998) proposed an obscu- 
ration threshold at Nznu = 1.4 x 10'-' cm"^ (above which the 
background QSO may be rendered invisible by dust obscura- 
tion), shown as a solid (blue) line in our Fig. 1. This was only an 
empirical limit based on the sample used by Boisse et al. (1998). 
It is thus not surprising that a few systems (with Zahs < 1-5) 
in Fig. 1 do lie above this threshold. Most of these systems 
were observed recently by Khare et al. (2004), Meiring et al. 
(2006a,2006b) and Peroux et al. (2006a, 2006b) and many of 
them have strong Mg II and/or Fe II lines in the SDSS spectra 
of the QSO's. Recently, Herbert-Fort et al. (2006) have reported 
the presence of two metal strong systems (with Zats > 1-6) in 
SDSS DR3, having Zn II column density greater than the obscu- 
ration threshold (SDSS1610-I-4724, log(Nz„ n)=13.4+0.03 and 
SDSS1709+3258, log(Nzn ii)=13. 19+0.03) and have estimated 
that ^ 5% of the SDSS-DR3 DLA population with Zabs > 2.2 in 
QSOs with r < 19.5 have similar Zn II column densities. 



2.2. Establishing the reality of the observed trend 

In this section we present some relevant observational facts and 
discuss the known selection effects to suggest that the Boisse 
plot may indeed not be affected by these and that the trend of 
decreasing abundance with increasing Nh i may be real. 

We first discuss the observed extinction in QSOs and point 
out that it is much smaller than that observed in typical Milky 
Way sight lines and may not cause significant dimming of the 
background QSOs. 

Vladilo & Peroux (2005) have shown that the extinction Ay 
is proportional to Nzn ii for Milky Way sight lines. Their relation 
is plotted in Fig. 2. Also plotted is a similar relation obtained for 
the sub-samples of Y06. It can be seen that a correlation does 
exist between Ay and Nzn ii for QSO absorbers but that the av- 
erage extinction per Zn ion in the these absorbers appears to be 
smaller than that in Milky way by a factor of three. 

The maximum value of E{B - V) found by Y06 for their sub- 
sample (#27) of most reddened systems was 0.085 while that for 
the sub-sample (#23) of systems with the rest frame equivalent 
width of Mg II A2196 < 2.0 A is as low as 0.003. The observer 
frame Ay, for the average redshift of 1.33 for their samples, as- 
suming a l/zl extinction law, is thus smaller than 0.6 and 0.02, 
respectively for the two samples. These values suggest that the 
dust obscuration by the absorbers observed towards the SDSS 
QSOs, may not be very important and may not cause significant 
decrease in the brightness of the QSOs. A similar conclusion is 
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Fig. 2. Absorber rest frame average extinction. Ay, assumed to 
be 2.93xE(B - V) vs. log(Nz„ii) for the Y06 results. The linear 
best fit relation Ay = O.lxlO'^ Nznii is shown as solid (red) 
line. The best fit line obtained for Milky Way sight lines. Ay 
= 0.3x10'^ Nznii by Vladilo and Peroux (2005) is shown, as 
dashed (black) line for comparison. 



drawn by Murphy & Liske (2004) from the analysis of spectral 
energy distributions of QSOs with DLAs at Zahs ~ 3 in the SDSS 
DRl. Ellison et al. (2005a) found no significant difference be- 
tween the B-K colours of radio selected QSOs with and without 
DLAs, indicating that dust obscuration is not very important. 

Below we discuss other indirect evidence indicating that dust 
obscuration is not significant in QSOs. 

Wolfe et al. (2005), from the observed depletion pattern of 
elements, estimated the fraction of dust obscured QSOs to be 
< 10%. Schaye (2001; also see Zwaan & Prochaska 2005) has 
shown that there is an upper limit of < 10~^ cm"^ on the neu- 
tral hydrogen column densities of DLAs because of conversion 
of H 1 to H2 and not because of dust obscuration. Herbert-Fort 
et al. (2006) have presented several arguments against the ob- 
scuration bias at redshifts > 1 .6, most importantly, one based on 
the magnitude distribution of the parent QSOs of metal strong 
systems. 

Samples 11 and 12 of Y06 had i magnitude smaller/larger 
than 19.1. Coincidently, this is the cutoff for the SDSS target- 
ing algorithm. Thus the bright sample consists of QSOs which 
should be unreddened to enter into the SDSS QSO catalog while 
the faint sample consists of candidates which were observed (and 
found to be QSOs) on the basis of their being X-ray sources, ra- 
dio sources etc. Thus one would expect the faint sample to be 
more reddened. However, both samples have similar values of 
E(B - V) and similar distribution of A(g - i) values. Thus there is 
no evidence for higher extinction in faint SDSS QSOs from the 
Y06 study. Prochaska et al. (2005) measured 40+20% higher 
gas density in DLAs towards bright QSOs than towards faint 
QSOs in SDSS DR3 having DLAs (z„;„ > 2). This is contrary to 
the obscuration bias and they suggest that gravitational lensing 
due to DLAs may be important. Vanden Berk et al. (1997) had 



found a significant excess of C IV systems in bright QSOs and 
interpreted this to be evidence for gravitational lensing. Some 
of these C IV systems could however be intrinsic to the QSOs 
(Richards 2001). Menard (2005) showed that lensing due to in- 
tervening Mg 11 systems with rest equivalent width smaller than 
1.5, could brighten QSOs by up to -0.2 magnitudes. Though the 
two eff'ects (lensing and dust obscuration) have diff'erent origins 
and need to be understood further, the results of Prochaska et al. 
(2005) indicate that the effects of lensing dominate over obscu- 
ration effects that may be present. 

The above arguments seem to indicate that dust obscuration 
is not important even in faint QSOs in the SDSS sample. We, 
however, note that all these arguments are based on samples of 
optically selected QSOs and can only apply if the dust content 
in QSO absorbers has a continuous distribution. We note that 
a few dusty DLAs have been observed (e.g. Junkkarinen et al. 
2004; Motta et al. 2002; Wild et al. 2006) but in none of these 
cases E{B - V) exceeds 0.42. Such values of E(B - V) may push 
bright QSOs below the cutoff (19.1) of SDSS targeting algorithm 
but these will still be present in the faint sample of Y06. Some 
empty fields have been observed in optical observations towards 
radio QSOs by Jorgenson et al. (2006) which could potentially 
be obscured by dusty absorbers, though the authors have argued 
against such a possibility on the basis of the DLA statistics ob- 
served in radio selected QSOs. 

Some theoretical studies (e.g. Vladilo and Peroux 2005) have 
argued that up to 30% to 50% of the QSOs may be missed 
in magnitude limited surveys due to dust obscuration. Smooth- 
particle-hydrodynamics simulations do need to introduce dust 
obscuration to explain the observed (low) DLA metallicity (Cen 
et al. 2003). Our arguments above do not rule out the possibility 
of a bimodal distribution of dust columns, with a population of 
dusty absorbers which have pushed the background QSOs below 
the observational limit of the optical surveys. Though this pos- 
sibility can only be verified with observations of fainter QSOs, 
no compelling evidence for such a population is found from the 
study of metallicity, column density distribution and mass den- 
sity of H 1 in DLAs in radio selected samples of QSOs (Akerman 
et al. 2005; Jorgenson et al. 2006). Their study suggests that the 
optically selected samples give a fair census of the population of 
DLA absorbers. It should, however, be noted that these samples 
are still small and cover redshifts > 1.86 only. 

From all the above arguments, it seems very likely that not 
many points in the upper right corner of the Boisse plot are 
missed due to dust obscuration. 

The missing points in the lower left corner of the Boisse plot 
have been attributed to an artifact of the sensitivity of observa- 
tions (Boisse et al. 1998). A few systems may indeed be present 
in this region as seen from the few upper limits there in Fig. 
1 . However, it may be noted that in the analysis of Y06, solar 
metallicity was estimated in the sub-sample (#24) of 698 sys- 
tems with average Nh i of about 10^" cm"^. Most (638) of these 
systems did not have detectable Zn 11 lines in their individual 
SDSS spectra. We note that even a solar metallicity system with 
Nh I of about 10^° cm"-, may not produce detectable Zn 11 lines 
in an individual SDSS spectrum which typically has S/N ~ 10- 
15 and a 3cr detection limit in the observer frame of > 0.3 A. So 
the absence of Zn 11 lines in an individual SDSS spectrum does 
not mean low Zn abundance. The composite spectrum of 698 
systems has S/N higher than that in a typical SDSS spectrum by 
a factor of ~ 25-30, enabling the detection of Zn 11 lines. In the 
sub-samples of Y06 (#s 21 and 22) comprised of systems with 
detectable Zn II and Cr 11 lines, the metallicity was estimated 
to be lower than that in the sample of 698 systems. It thus ap- 
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pears that their metallicity measurements, being averages over 
large samples, are not affected by observational htnitations in 
detecting weak lines. We also note that the average solar metal- 
licity in the sub-sample # 24 can not be an artifact of a possible 
high metallicity tail of the metallicity distribution in sub-DLAs. 
The composite spectrum, used for abundance determination is 
the arithmetic mean of individual spectra. Thus, if 90% of the 
698 systems had subsolar abundance with [Zn/H]=-1 and only 
10% had solar abundance then, assuming an H 1 column density 
of 10^° cm"^ for the systems, the average equivalent width of 
A2026 line of Zn II will be 0.9 x 0.007 -i- 0.1 x 0.07 = 0.013A, 
which will correspond to highly sub-solar abundance of Zn 11 
for this sub-sample as opposed to the solar abundance derived 
by Y06. It may also be noted that the sub-samples #23 and 24, 
do have a large number of systems with low Nh i which in princi- 
ple, could have had low [Zn/H] and could then have been present 
in the lower left corner of the Boisse plot. Sub-sample #23 which 
is the sample of systems with Wjvig u < 2.0 A, includes all the 
systems in sub-samples #s 2, 3 and 4 (with WMg n < 1.91 A, all 
together 397 systems) which have E(B - V) < 0.001 and thus 
Nh I < 5 X 10^'^ cm"^. Similarly sub-ample #24 has 698 systems. 
Even assuming it has all the systems with Wjvig u > 1.91 A (i.e. 
all the systems not included in sub-amples #s 2, 3 and 4) from 
the original sample of 809 systems, it will still have > 286 sys- 
tems with Nh I < 5 X lO'^ cm"^. However, the high values of the 
average abundances for these sub-samples suggest that most of 
these low Nh i systems do not have low abundance and thus do 
not occupy the lower left corner of the Boisse plot. We also note 
that it is possible that these two sub-samples do include some 
DLAs. However, as argued above, the DLAs do not have high 
abundances and can only reduce the average abundances of the 
sub-samples, thereby, making the case for higher abundance of 
sub-DLAs even stronger. 

On the conservative side, we point out that even if some 
points (in TLS) are indeed missing in the lower left region of 
Fig. 1 due to sensitivity limit of observations, the average metal- 
licity at low N H I (in TLS) will still be higher than that at high 
N H I due to the large number of high metallicity, low N h i sys- 
tems that have been observed (assuming the obscuration bias to 
be absent, as argued above). It thus appears that the sub-DLA 
metallicity is indeed higher than that of DLAs and is close to the 
solar value, at low redshifts. 

It thus seems very hkely, that the observed trend of decreas- 
ing abundance with increasing Nh i is not due to selection effects 
and is real and that DLAs are not among the major metal carri- 
ers in the universe. We note that super-solar abundances have 
been observed in seven sub-DLAs (Pettini et al. 2000; Khare et 
al. 2004; Prochaska et al. 2006; Meiring et al. 2006b). Near so- 
lar abundance has been estimated in one NDLLS (Jenkins et al. 
2005). Bergeron et al. (1994), estimated abundances ([X/H]), of 
several Mg 11 Lyman Limit systems with low Lyman limit opti- 
cal depth, with redshifts between 0.1 and 1.1, to be between -0.5 
and -0.3 dex. Super-solar abundances have also been estimated 
for three systems with redshifts between 0.7 and 0.95 (Charlton 
et al. 2003, Ding et al. 2003; Meseiro et al. 2005) and near-solar 
abundance has been estimated in one system at redshift of 0.064 
(Aracil et al. 2006); aU these systems have Nh i < 10^'' cm~^. We 
however caution that the derived abundance in NDLLS and sys- 
tems with smaller H 1 column densities do depend on the details 
of photoionization calculations. 

Having argued above, for the reahty of the observed depen- 
dence of Zn metallicity on H 1 column density, below, we con- 
sider its implications. In particular, we consider consequences 



of the hypothesis that the sub-DLAs, and possibly NDLLS, have 
higher abundances than DLAs and represent a selection of galax- 
ies that are the major metal carriers in the universe. 

3. Implications 

Recently a mass-metallicity relation has been discovered by sev- 
eral groups. Tremonti et al. (2004), from the imaging and spec- 
troscopy of 53,000 star-forming galaxies at z ~ 0.1, found a 
tight (±0. 1 dex) correlation between stellar mass and metallic- 
ity spanning over 3 orders of magnitude in stellar mass and a 
factor of 10 in metallicity. Savaglio et al (2005), from a sam- 
ple of 56 galaxies identified a strong correlation between mass 
and metallicity at 0.4 < z < 1.0. They predict that the gener- 
ally metal poor DLA galaxies have stellar masses of the order of 
10^-^ Mo (with a dispersion of 0.7 dex) from z=0.2 to z=4. Erb et 
al. (2006) have obtained a mass-metallicity relation in a sample 
of 87 galaxies at < z > ~ 2.26 which is similar to the relation 
for local galaxies (Tremonti et al. 2004) except for an offset of 
0.3 dex in metalhcities, indicating that galaxies of a given mass 
have lower metallicity at high redshift. They however, note that 
the uncertainty in the metallicity offset between the z ~ 2 and lo- 
cal galaxies is approximately a factor of 2, about the same as the 
offset itself. Ledoux et al. (2006) found a correlation between the 
metallicity and velocity widths of lines of low ionization species 
over two orders of magnitude in metallicity, at 1.7 < z < 4.3. 
Assuming velocity widths to be a measure of mass, their mass- 
metallicity relation is consistent with that found for local galax- 
ies by Tremonti et al. (2004). We, however, note that Bouche et 
al. (2006) found an anticorrelation between halo mass and Mg II 
A 2796 equivalent widths which seems to go against the results 
of Ledoux etal. (2006). 

At low redshifts the metallicity of most DLAs is almost an 
order of magnitude lower than the solar value, while, as sug- 
gested by results of Y06 and those of Kulkarni et al. (2006b), the 
mean sub-DLA metallicity seems to be close to the solar value. 
The mass-metallicity relation of Tremonti et al. (2004) would im- 
ply stellar masses of about 10^^ Mq and < 10^ M© respectively 
for the sub-DLAs and DLAs. These numbers will be smaller by a 
factor of 2 if the mass-metallicity relation of Erb et al. (2006) 
is used. If the stellar metallicity is lower/higher than the gas 
phase metallicity (as determined by the abundances in DLAs and 
sub-DLAs) the masses of the DLA and sub-DLA galaxies may 
be correspondingly lower/higher, however, the ratio of the two 
masses will not be affected. 

We thus propose that the sub-DLAs and DLAs in general 
represent different types of galaxies. By and large, the sub-DLAs 
are produced by massive galaxies, with higher metalhcity, while 
the DLAs are produced by less massive galaxies with lower 
metallicity. A few sub-DLAs may indeed arise from lines of 
sight which encounter low H I columns through DLA galax- 
ies and therefore, have lower metallicity. These will give rise 
to points in the lower left region of Fig. 1. Similarly, a few 
DLAs may arise in lines of sight through sub-DLA galaxies and 
may have higher metalhcities though the probability of this hap- 
pening may be small. In our Galaxy, for instance, clouds with 
log(NH i)>20.3 are very small and represent only a tiny part of 
the cross section of the entire Galaxy. 

It has been suggested that the low metallicity found in DLAs 
is caused by metallicity gradients. Differences between emis- 
sion fine and absorption line metaUicities have been observed 
in a few cases (Chen et al. 2005; Ellison et al. 2005b) but not 
in others (Schulte-Ladbeck et al. 2005; Bowen et al. 2005). The 
metallicity gradients in nearby spirals are fairly weak (Bresolin 
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et al. 2004) and can not explain the low metallicity in DLAs 
with small impact parameters (e.g. Kulkarni et al. 2005). It has 
been suggested in several studies (e.g. Kauffmann 1996; Das et 
al. 2001; Zwaan et al. 2005) that DLAs may result for small val- 
ues of impact parameters while larger impact parameters through 
the same absorbers may give rise to NDLLS and Lyman a for- 
est systems. In this case, abundance gradients can not be in- 
voked to explain lower abundances in DLAs as compared to 
sub-DLAs/NDLLS. This scenario has also not been verified ob- 
servationally. 

Kulkarni et al. (2006b) have shown that the Nh i weighted 
mean metalUcity for sub-DLAs is a factor of 6 higher than that 
for DLAs at low redshifts. They estimate that at these redshifts, 
the contribution of the ISM in sub-DLAs to the cosmic metal 
budget may be several times that of ISM in DLAs. Prochaska et 
al. (2006), based on their observations of two metal strong sub- 
DLAs (Super LLS in their paper) and results from their sub-DLA 
survey, estimated that the ISM in sub-DLAs may contribute at 
least 15% to the metal budget of the universe at z ^ 2. 

If the galaxies that give rise to sub-DLAs are indeed more 
massive than the galaxies that give rise to DLAs, as suggested 
above, then an individual sub-DLA galaxy will contribute more 
to the stellar mass of the universe than an individual DLA 
galaxy. Also, in CDM cosmology, as the massive galaxies form 
by mergers of smaller galaxies which triggers star formation 
leading to a higher rate of metal production, they are not only 
expected to show a higher rate of metallicity evolution (as is 
found by Kulkarni et al. 2006b) but are also expected to con- 
tribute more mass at lower redshifts. It is thus possible that the 
contribution of sub-DLA producing galaxies to the cosmic metal 
budget at lower redshifts may indeed be considerably higher 
than 15% (as estimated by Prochaska et al. (2005) at z — 2) 
and may help alleviate the missing metals problem. 

Boissier et al. (2003) compared the observed properties of 
DLAs with the predictions of simple models of present day disk 
galaxies and showed that low surface brightness galaxies con- 
tribute as much as spirals to the number and H 1 mass of DLAs. 
Zwaan et al. (2005) have shown that properties of DLAs are con- 
sistent with their forming in galaxies of various morphological 
types, with 87% of DLA cross-section being contributed by sub- 
L, galaxies. Semi-analytic models indicate sub-L, galaxies to be 
major contributors to DLA cross-section (Okoshi & Nagashima 
2005). Thus there seems to be some evidence for a significant 
fraction of DLAs to be associated with sub-L, galaxies. 

Deep imaging has shown low redshift DLAs to be often as- 
sociated with dwarf galaxies (Rao et al. 2003). Le Brun et al. 
(1997) suggested that the DLA galaxies (which were selected 
on the basis of a damped Lyman a line or 21 cm absorption 
or a very high Mg II/Fe II ratio) strongly differ from the Mg 
11 selected galaxies, the latter being mostly luminous galaxies 
with evidence of recent star formation activity. Indeed emission- 
line imaging searches suggest that a large fraction of DLAs ap- 
pear to have low star formation rates (Kulkarni et al. 2006a and 
references therein). Deep imaging of three sub-DLA galaxies at 
z < 0.7 shows them to be associated with L> 0.6 L,, disk/spiral 
galaxies (Zwaan et al. 2005). Note that the Unes of sight through 
our Galaxy are, statistically speaking, mostly metal rich sub- 
DLAs with Nhi < 2 X 10^" cm"^, perpendicular to the plane. 
All these support our hypothesis about the nature of DLA and 
sub-DLA galaxies. Chen and Lanzetta (2003), however, found 
the DLA galaxies to have a variety of morphologies. Djorgovski 
et al. (1996) and MoUer et al. (2002, 2004) have observed a few 
high-redshift (z >2) absorbing galaxies in Lyman a emission and 
found their properties to be similar to disk/Lyman break galax- 



ies. Chun et al. (2006), using adaptive optics imaging for 4 DLAs 
and 2 sub-DLAs at z < 0.5, found most of the candidate absorber 
galaxies or their companions to have low-luminosity ( < 0.1 L.). 

Thus, it seems that the results of some of the DLA and 
sub-DLA imaging studies are consistent with the hypothesis we 
make here. More systematic imaging surveys are needed to con- 
firm the ideas presented in this paper. Deep imaging in K-band 
for absorbers with high abundances at Zats < 1 .5 should be able 
to confirm the presence of massive (red) galaxies, while imaging 
searches in narrow, optical emission lines should help in con- 
firming that DLAs are mostly dwarf galaxies. 

4. Conclusions 

We have studied the dependence of metalUcity on Nh i in QSO 
absorbers with Nhi > lO" cm"-^ and have discussed various 
selection effects that may be giving rise to the observed trend. 
We have argued that the selection effects are not important and 
that the observed trend is real. Our conclusions based on these 
arguments, and subject to confirmation by future observations, 
are as follows: 

1 . The amount of dust in QSO absorbers is small and is not re- 
sponsible for missing many QSOs in magnitude limited sur- 
veys. We can not, however, rule out the possibility of a bi- 
modal distribution of dust columns such that there may exist 
a population of dusty absorbers which push the background 
QSOs below the observational threshold of current optical 
spectroscopic studies and is completely invisible. 

2. The metallicity in QSO absorbers with Nhi > 10'^ cm'^, 
decreases with increase in H I column density of these ab- 
sorbers. The trend possibly continues to lower H I column 
densities. The sub-DLAs thus have higher metal abundances 
as compared to the DLAs at redshifts between and 2. 

3. The observed mass-metallicity relation suggests that most 
DLAs are associated with low mass (< lO' Mq) galax- 
ies while most sub-DLAs are associated with massive spi- 
ral/elliptical galaxies. It is possible that the non-DLA LLS 
may also be metal rich and may be associated with massive 
galaxies. 

4. The sub-DLA galaxies will contribute a larger fraction of to- 
tal mass (stellar and ISM) and therefore metals, to the cosmic 
budget, specially at low redshifts, as compared to the DLAs. 
The Sub-DLAs and possibly, non-DLA LLS galaxies, may 
contain a much larger fraction of the metals at z < 1 than has 
been appreciated. 

5. The few imaging studies, of galaxies responsible for quasar 
absorption fine systems done so far, are ambiguous on the 
morphology of DLA and sub-DLA galaxies. More system- 
atic, deep imaging in r-band, Ks-band and in narrow emis- 
sion lines is essential to confirming the inferences of this pa- 
per. 
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